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We study the effects of a magnetic impurity on the behavior of a S = 1/2 spin chain. At T = 0, 
both with and without an applied uniform magnetic field, an oscillating magnetization appears, 
whose decay with the distance from the impurity is ruled by a power law. As a consequence, 
pairwise entanglement is either enhanced or quenched, depending on the distance of the spin pair 
with respect to the impurity and on the values of the magnetic field and the intensity of the impurity 
itself. This leads us to suggest that acting on such control parameters, an adiabatic manipulation of 
the entanglement distribution can be performed. The robustness of our results against temperature 
is checked, and suggestions about possible experimental applications are put forward. 
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The prospect of quantum technology, i.e., the design 
and realization of quantum devices, is triggering an in- 
creasing interest in the manipulation of quantum systems 
in view of several applications, amongst which is quan- 
tum information processing [ij. Despite the large vari- 
ety of experimental proposals, the relevant properties of 
most of such devices can be adequately described by a few 
classes of models. Quantum low-dimensional spin mod- 
els constitute one such class: They may describe different 
physical systems and, at the same time, their Hamilto- 
nian is simple enough to allow an easy recognition of the 
parameters responsible for specific effects, so that simple 
controls may be designed to drive focused manipulations. 
Sticking to strictly one-dimensional models, a possibility 
of manipulation is offered by the introduction of local im- 
purities, either magnetic or not, which break the transla- 
tional symmetry. In the last decade, several studies have 
been devoted to low-dimensional s^tems with inhomo- 
geneities 0, i i, S, i, 0, i, & Ell, [II, revealing 
peculiar phenomena which occur in the vicinity of the 
impurity. 

In this paper, we consider the less investigated case of 
magnetic impurities, which can couple both with neigh- 
boring spins and/or with external magnetic fields, offer- 
ing a remarkable control flexibility. In particular, aimed 
at possible application in the realm of quantum compu- 
tation, we consider a single impurity in the 3 = 1/2 XX 
chain in a transverse field in its quasi-long-range ordered 
phase. In the absence of impurity and due to the high 
symmetry of the exchange interaction, such model dis- 
plays values of bipartite entanglement larger than those 
observed in less symmetric models 13|. The impurity is 
modelled in terms of an additional field located at one 
precise site of the chain. We study the ground-state 
properties and thermodynamic behavior when the field 
and/or the magnetic strength of the impurity are varied. 



At zero temperature, we find an infinite penetration 
depth for the effects of the magnetic impurity, namely 
it produces spatial modulations of the spin configuration 
which decay polynomially with the distance from the im- 
purity itself. A staggered magnetization and an oscillat- 
ing susceptibility result, as a most general consequence of 
translational symmetry breaking [2]. One striking effect 
of the emergence of such spin structure is the local en- 
hancement of pairwise entanglement: When no uniform 
field is applied, the concurrence [13, m of the spm-pair 
next to the impurity doubles with respect to the transla- 
tional invariant case and, consequently, an entanglement 
re-distribution takes place all over the chain. With finite 
magnetic field, there always exists a most entangled spin 
pair, which can be moved along the chain by tuning ei- 
ther the field or the impurity strength. These parameters 
can thus play the role of "knobs" for the manipulation 
of both the magnetic and the entanglement properties of 
the chain. Furthermore, as they appear in terms that 
commute with the Hamiltonian, adiabaticity is guaran- 
teed while tuning their values. The essential features of 
the above picture, and in particular the possibility of se- 
lectively enhance pairwise entanglement, survive at finite 
temperature. This gives a concrete possibility to experi- 
mentally test our findings. 

Let us consider the Hamiltonian 



n 

J 



E 



-a a. 



, 



(1) 



where cr^ (e = x, y, z) are the the Pauli operators for 
the spin at site n, and periodic boundary conditions 
are assumed. A uniform magnetic field H is applied 
along the quantization axis, which takes values such that 
h = g^ji^H/J e [—1,1]; the exchange integral J is set 
to unity throughout the paper. The magnetic impurity 
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is located at site and represented by a local magnetic 
field a [31 . As Ti. can describe several physical systems 
beyond the mere spins (e.g. atoms loaded in optical lat- 
tices), one may think of correspondingly several way of 
introducing a localized defect. In general the impurity 
breaks the translational symmetry and as a ^ oo it ren- 
ders the system equivalent to an open-end chain. 

The Hamiltonian ([1]) can be diagonalized by standard 



methods [17|, |18|, resulting in an energy band and in a 



localized level 
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Here, = h — (sgna)-\/l -I- is the energy of the lo- 
calized state, while the corresponding fermion operator 
is a combination of the Jordan- Wigner ones {cj^} 
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with ^ ^ = — ln(Vl — — \a\). Concerning the band, 
Ek=h—cos9k is the unperturbed dispersion relation 
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The impurity is seen to produce i) the appearance of 
a discrete energy level (lying above or below the band 
depending on the sign of a), exponentially localized in 
space with a characteristic length ^; and ii) the distor- 
tion of the energy-band states, whose amplitudes, besides 
the plane wave term, contain a back-scattering contribu- 
tion due to the defect. It is precisely the interference 
between these two terms that gives rise to the oscillating 
patterns discussed below in both the magnetization and 
the susceptibility. 

Using the Wick's theorem, we calculate one- and two- 
points correlation functions which, due to the lack of 
translational symmetry, cannot be expressed as Tocplitz 
determinants, so that their expressions become more and 
more involved with the relative site distance. We find: 

(a,X^) = l-2^(|v..,,| 

(4) 



where is the set of filled states. Rotational symmetry 
gives (ajf) = (a-f) ^ and {a^a^^^) = {a^a^^g)] more- 
over, reflection symmetry with respect to site allows us 
to consider only r > in the following. 
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FIG. 1: (cr^> at ft = (left) and h = 0.5 (right), for a = 
(continuous line), a = 1 (filled square) and a — 5 (open 
circles). The insets show the local susceptibility for a — 5. 



In the thermodynamic limit, the local magnetization 
at zero temperature takes the form 
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with 6p = cos ^ h and where (cr^) = 1 — gives the 
result without defect, {a^)x = ~2\ipr,x\'^ is due to the lo- 
calized state (this term is present only for a > 0), while 
the last term is due to the interference mechanism men- 
tioned above. The essential finding embodied in Eq. ([5]) 
is the emergence of an oscillating magnetization in the 
vicinity of the impurity, as shown in Fig. [H 



Several works 0, i, i, i, @, 0, S, S,]lfl [Hi have ev- 
idenced that finite staggered moments are induced by 
open ends in antiferromagnetic spin chains and that the 
spatial decay of the related magnetization is ruled by 
the (bulk) correlation length of the infinite, translational 
invariant, system. Previous analysis, however, referred 
to spin chains in disordered phases, with exponentially 
decaying correlations, cither because of the finite tem- 
perature or due the symmetry of the Hamiltonian itself, 
as in the case of the Haldane chains. On the contrary, 
our model is independent of the sign of the exchange 
coupling and; at T = 0, it is in a quasi-long-range or- 
dered phase with an infinite correlation length and a 
much slower (power law) decay of the correlation func- 
tions: The effects of the impurity are consequently felt 
by more distant spins along the chain. Indeed, there ex- 
ists a precise analogy Q between our model and the case 
of a one-dimensional fermionic system in the presence of 
one impurity, where Friedel oscillations occur [l9| : For 
a ^ 1 and r ^ 1 Eq. ([5]) gives the same ^ behavior 
predicted for the oscillations of the fermion density in a 
Luttinger liquid with a defect [23|. It is noteworthy that, 
at variance with the case of the XX chain with locally 
inhomogeneous exchange interaction 0], the presence of 
the magnetic impurity gives rise to Friedel-like oscilla- 
tions even when no uniform field is present. 

In general, for a fixed h, the value of a determines the 
absolute value of the local magnetization, while, indepen- 
dently of a, the field fixes the spatial periodicity p of the 
oscillations (sec Fig. [1]). In fact, we see that p coincides 
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with the period of the connected correlation functions 
of the translational-invariant model, i.e. P = [HI- 
This behavior is understood by observing that the local 
switching-on of (tr^) under the effect of the local field at 
site is mediated by the generalized magnetic suscepti- 
bility, which in turns has the same spatial dependence of 
the correlation functions. 

Let us now consider the local magnetic susceptibil- 
ity. From Eq. (P, we get Xr{a,h) = 2/{tt^/1 - K^) + 
Xr'*(a, h), where the alternating term is 
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FIG. 2: Concurrence C12 (left) and C23 (right) vs h, for sev- 
eral values of a. 



The superposition of a uniform term and a spatially oscil- 
lating one, is fully analogous with the results relative to 
the spin-i Hcisenberg antiferromagnctic chain with open 
ends or bond impurities 0. The alternating term 
survives throughout the chain with a spatial structure 
that, being r an integer, follows peculiar, bcating-likc, 
patterns, depending on the value of h. When h ~ 0, one 
has Xr(ck, 0) = § ^1 — ^~i]^Ji ^ meaning that, for a » 1, 
odd-labelled spin increase their susceptibility up to the 
value while even-ones have zero susceptibility. 

Noticeably, Xr{o:,h) does not diverges at the critical 
point h = I. In fact, while in the translational invari- 
ant model the susceptibility diverges as (1 — h)~^^'^ for 
/i — !■ 1, wc find here that in the same limit Xr{ct,h) ~ 



2^2(1-/1) / 1+2™ 



As a consequence of the transla- 



1^+2(1-/1)^ 

tional symmetry breaking, singularity is cancelled and all 
susceptibilities vanish at the critical point. 

The spatial modulation of the magnetic susceptibility 
described in Eq. ^ , offers a chance to selectively act on 
local spins. This possibility is at the hearth of the fol- 
lowing discussion about local enhancement and possible 
transfer of entanglement. 

As stated in the introduction, due to the high sym- 
metry of its Hamiltonian, the XX model has a markedly 
entangled ground state, with next neighbor concurrences 
larger than those observed in less symmetric models [l^ . 
Furthermore, the rotational invariance forces bipartite 
entanglement to be of antiparallel type [13, H^. This 
property holds also when translational invariance is bro- 
ken. This fact, together with the emergence of the stag- 
gered structure described above, suggests an enhance- 
ment of the nearest neighbor concurrences. This is shown 
to occur in Fig. [21 where C12 and C23 are displayed vs 
h, for different values of a. In particular, at zero field 
and for large a, C12 is more than twice its value without 
defect, and a similar enhancement is observed in C23 for 
two different values of h. In general, Cr,r+i displays r 
maxima, whose position in the [—1,1] field interval is in- 
creasingly symmetric, with respect to the value /i = 0, as 
a increases. The maximum values of Cr,r+i are always 
larger than those attained for a = 0, but the difference is 



smeared out with increasing r (i.e., by moving far from 
the impurity). The highly non trivial behavior of Cr^r+i 
with respect to h, r and a, is formally due to the inter- 
ference phenomenon embodied in Eq. ([3]). Physically it 
is a consequence of an entanglement re-distribution, re- 
sulting from the constraints imposed by the monogamy 
inequality [2J, [25[ on a locally modulated spin-structure. 

In particular, for h = and independently of a, the 
most entangled spin pair is the one adjacent to the im- 
purity and the one-tangle = 1 — (c^)^, with r 7^ 0, 
does not vary much with r and it is substantially insensi- 
tive to a. The observed alternating behaviour of C^^r+ij 
which gets a larger or smaller value with respect to the 
translational invariant case according with r being odd 
or even, respectively, can thus be interpreted as a con- 
sequence of the monogamy inequality. Furthermore, the 
difference between Cr,r+i and Cr+i,r+2 increases with a 
and a sort of entanglement dimerization is observed for 
large impurity strength. 

For h ^ 0, Cr.r+i displays an increasingly complicated 
dependence on h and a, which gives rise to a sharp- 
cut mechanism of entanglement transfer for /i < 1. In 
such case indeed, as shown in Fig. [3] for h = 0.9 and 
0.99, a most entangled spin pair is still detected at a dis- 
tance from the impurity that increases with h. Hence, 
by varying either the defect strength or the magnetic 
field, one can drive the maximum of the concurrence 
along the chain, meanwhile reducing the entanglement 
between spin pairs next to the defect. It is of absolute 
relevance that the couplings of the spin chain with both 
the uniform and the local fields are described by terms 
which commute with the bare exchange XX Hamiltonian. 
This means that a time variation of the control parame- 
ters h and a induces a fully adiabatic dynamics as non- 
equilibrium configurations are never accessed. 

When temperature T is switched on, thermodynamic 
averages are given by Eqs. ([4]) weighted by fermionic pop- 
ulation factors. As seen in Fig.[31 the spatial distribution 
of the magnetization is quite robust against temperature, 
essentially because of the energy gap between the local- 
ized level and the bottom of the band. This is also un- 
derstood by noticing that the magnetic moment at the 
impurity, (ctq), keeps finite at finite temperature. In fact, 
neglecting the contribution of the band amplitudes, we 
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FIG. 3: Concurrence Cr,r+i vs r at = 0.9 (left), h = 0.99 
(right) for several values of a. 
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FIG. 4: (a^) vs T for h = and a = 1 (left); concurrence 
Cr,r+i at h — 0, a = 5 VS T (right). 
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■■ tanh 



where the functional 



'0/ - VT+^ """" 2T 
dependence of the hyperbolic tangent on T accounts for 

the temperature robustness. 

As for the concurrences, a marked dependence on r 
survives at T ^ and thermal entanglement is found for 
selected spin pairs, with C^.r+i displaying a maximum at 
finite temperature for even r, while it decreases monoton- 
ically for r odd. Although several systems which develop 
thermal entanglement are known, our model provides the 
additional possibility to select the spin pair for the en- 
hancement (reduction) of entanglement. 

Possible experimental studies of the above issues set in 
several frameworks. First of all, quasi-one dimensional 
magnetic compounds properly doped so as to introduce 
diluted magnetic impurities, might be considered [2^ (to 
this respect, we verified that our results still hold in the 
case of two impurities, provided they are far enough from 
each other) . In this case NMR imaging of the local mag- 
netization near impurities could be performed, in analogy 
with the case of non- magnetic impurity the same 

technique couldbe used to probe the spatially modulated 
susceptibility 
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Moreover, muon-spin relaxation ex- 
periments have been extensively performed to investi gate 
the effects of magnetic doping in high-Tc cuprates [28[ 
and might also be proper tools for testing our findings, 
since the muon acts as both the impurity and the probe 
for the local magnetic arrangement [2^. In this frame- 
work, however, no external tuning of the parameter a 
is possible, as its value would result from the magnetic 
moment of the impurity, uniquely fixed by the dopant 
properties. Nonetheless, control upon the uniform field 
would still be feasible. Another possible implementation 
is suggested by the recent proposal js^l for realizing effec- 



tive spin systems by coupled micro-cavities: in this case, 
the uniform magnetic field might be tuned by varying the 
frequencies of the driving lasers, while the impurity term 
could result from the level structure of the micro-cavity 
representing the spin sitting at the site of the impurity. 
Atomic condensates in optical lattices can also be consid- 
ered, as it has been shown that they can simulate bosonic 
and fermionic models which can be mapped to the XX 
model in a transverse field in its quasi-long-range ordered 
phase [iH . 
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